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ABSTRACT 

We  present  a  detailed  multi-epoch  analysis  of  31  potential  Southern  hemisphere  radio  cal¬ 
ibrators  that  were  originally  observed  as  part  of  a  programme  to  maintain  the  International 
Celestial  Reference  Frame.  At  radio  wavelengths,  the  primary  calibrators  are  active  galactic 
nuclei  (AGNs),  powerful  radio  emitters  which  exist  at  the  centre  of  most  galaxies.  These 
are  known  to  vary  at  all  wavelengths  at  which  they  have  been  observed.  By  determining  the 
amount  of  radio  source  structure  and  variability  of  these  AGNs,  we  determine  their  suitability 
as  phase  calibrators  for  long  baseline  radio  interferometry  at  2.3  GHz.  For  this  purpose,  we 
have  used  a  set  of  complementary  metrics  to  classify  these  31  southern  sources  into  five 
categories  pertaining  to  their  suitability  as  very  long  baseline  interferometry  calibrators.  We 
find  that  all  of  the  sources  in  our  sample  would  be  good  interferometric  calibrators,  and  almost 
ninety  per  cent  of  them  would  be  very  good  calibrators. 

Key  words:  techniques:  high  angular  resolution  -  techniques:  interferometric  -  galaxies: 
active  -  galaxies:  jets  -  quasars:  general  -  radio  continuum:  galaxies. 


1  INTRODUCTION 

High  angular  resolution  observations  of  weak  radio  sources  (where 
self-calibration  is  not  possible)  require  calibrator  sources  for  cor¬ 
rection  of  systematic  effects  and  effects  of  the  atmosphere  on  the 
measured  visibilities.  Atmospheric  fluctuations  cause  perturbations 
in  visibility  phase  which,  if  not  corrected,  seriously  limit  both  the 
sensitivity  and  image  quality  of  an  interferometric  array.  Phase  cal¬ 
ibrators  are  also  required  for  astrometric  observations.  An  ideal 
calibrator  would  look  the  same  on  all  observing  baselines.  It  should 
be  bright,  unresolved  or  at  least  compact  and  should  not  vary.  A 
calibrator  source  should  also  be  separated  from  the  target  sources 
by  as  small  an  angle  as  possible  in  order  to  look  along  the  same 
line  of  sight  throughout  the  atmosphere.  Therefore,  it  is  desirable  to 
have  calibrator  sources  evenly  distributed  across  the  whole  sky.  In 
practice,  at  radio  wavelengths,  calibrators  are  mostly  active  galactic 
nuclei  (AGNs),  whose  fundamental  source  of  power  is  believed  to 
be  the  accretion  of  matter  on  to  a  super-massive  black  hole  (e.g. 
Rees  1997).  They  are  known  to  vary  at  every  wavelength  at  which 
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they  have  been  studied.  AGNs  are  very  compact  and  isotropically 
distributed  around  the  sky.  They  are  also  very  distant  objects  and 
therefore,  generally,  have  no  discernible  proper  motions  on  the  sky. 
It  is  these  qualities  that  make  them  suitable  as  calibrators  at  radio 
wavelengths. 

Due  to  the  limited  number  of  radio  telescopes,  few  surveys  for 
calibrators  have  been  carried  out  in  the  Southern  hemisphere  (Fey 
et  al.  2004a, b,  2006;  Ojha  et  al.  2004a,  2005).  Hence,  there  are 
fewer  known  calibrators  in  the  south.  A  major  expansion  of  ra¬ 
dio  astronomy  observing  capability  is  underway  in  the  South¬ 
ern  hemisphere.  Two  Square  Kilometre  Array  (  SKA)  precur¬ 
sors,  the  South  African  Karoo  Array  Telescope  (MeerKAT;  Booth 
et  al.  2009)  and  Australian  SKA  Pathfinder  (ASKAP;  Johnston 
et  al.  2008)  are  presently  under  construction,  leading  to  the  SKA 
itself.  It  is  clear  that  interferometry  and  very  long  baseline  interfer¬ 
ometry  (VLB I)  in  the  Southern  hemisphere  need  a  dense  network 
of  calibration  sources  at  high  resolutions  and  a  range  of  frequen¬ 
cies.  MeerKAT  is  the  South  African  SKA  demonstrator  telescope. 
When  completed,  it  will  be  the  largest  radio  telescope  array  in 
the  Southern  hemisphere.  It  will  operate  in  two  frequency  ranges: 
0.58-2.5  and  8-14.5  GHZ  (Booth  et  al.  2009).  MeerKAT  will  par¬ 
ticipate  in  YLBI  observations  with  the  European  VLBI  network 
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and  other  VLBI  arrays.  The  SKA  will  operate  in  the  frequency 
range  0.7-25  GHz  (Carilli  &  Rawlings  2004;  Schilizzi  2007)  and  is 
expected  to  have  much  longer  baselines  for  which  a  southern  VLBI 
calibrator  list  will  be  essential. 

The  United  States  Naval  Observatory  (USNO)  in  collaboration 
with  National  Aeronautics  and  Space  Administration  (NASA),  the 
Laboratoire  d’Astrophysique  de  Bordeaux  (LAB)  and  the  National 
Radio  Astronomy  Observatory  (NRAO)1  has  since  1994  been  ob¬ 
serving  AGNs,  once  every  two  months,  through  the  research  and 
development  of  VLBI  (RDV)  experiments.  This  has  led  to  a  wealth 
of  data  of  sources  observed  simultaneously  at  2.3  and  8.4  GHz  at 
regular  epochs.  Images  of  these  AGNs  form  the  USNO’s  Radio 
Reference  Frame  Image  Database  (RRFID)  and  LAB’s  Bordeaux 
VLBI  Image  Database  (BVID)2  (Fey,  Clegg  &  Fomalont  1996;  Fey 
&  Chariot  1997,  2000;  Collioud  &  Chariot  2009).  The  RRFID  is 
a  data  base  of  about  6700  images  of  over  700  AGN  sources  com¬ 
piled  from  geodetic  and  astrometric  VLBI  experiments.  The  BVID 
contains  over  1800  images  of  824  radio  sources.  It  is  the  goal  of 
the  ongoing  RDV  programme  to  image  the  radio  reference  frame 
sources  on  a  regular  basis  to  monitor  them  for  variability  and  struc¬ 
tural  change.  Piner  et  al.  (2007)  discuss  jet  kinematics  in  a  subset 
of  the  RRFID  sources,  concentrating  on  the  8.4  GHz  observations. 
Their  survey  is  made  up  of  RRFID  sources  that  have  been  observed 
at  three  or  more  epochs  from  1994  July  to  1998  December. 

In  this  paper,  we  focus  on  2.3  GHz  observations  of  a  sample  of 
31  potential  southern  VLBI  calibrators  (with  declinations  between 
0°  and  —60°)  from  the  kinematic  survey  of  Piner  et  al.  (2007).  The 
choice  of  2.3  GHz  was  made  because  both  MeerKAT  and  ASKAP 
will  be  operating  at  frequencies  close  to  this.  We  seek  to  charac¬ 
terize  this  selection  of  southern  radio  sources  and  determine  their 
suitability  as  calibrators  for  southern  VLBI  experiments,  especially 
those  using  MeerKAT  and  the  SKA  when  completed. 

In  the  following  sections,  we  will  describe  the  observations  and 
data  reduction  process.  We  will  then  describe  the  imaging  and 
model-fitting  process  and  go  on  to  analyse  the  data  and  discuss 
the  results. 


2  OBSERVATIONS  AND  DATA  REDUCTION 

RDV  experiments  are  carried  out  using  the  10  antennas  of  the 
NRAO’s  Very  Long  Baseline  Array  (VLB A)  and  up  to  10  other 
antennas  across  the  globe,  including  Hartebeesthoek  (South  Africa) 
when  available.  The  use  of  the  global  array  greatly  improves  the 
u  —  v  coverage.  Hartebeesthoek  greatly  improves  the  u  —  v  coverage 
for  sources  in  the  south,  which  make  up  the  sample  in  this  paper. 

The  31  southern  sources  included  in  our  sample  are  listed  in 
Table  1  along  with  their  optical  properties.  The  source  sky  distribu¬ 
tion  is  shown  in  Fig.  1.  The  RDV  observing  epochs  are  from  1994 
July  to  2008  January,  giving  a  total  of  32  epochs  and  an  average 
of  20  epochs  per  source.  About  100  sources  were  observed  in  each 
24  h  observing  run,  with  an  average  on-source  time  of  15  min.  The 
on-source  time  is  not  continuous,  but  divided  into  scans  of  between 
one  and  several  minutes  long,  spaced  in  time  to  give  optimal  u  —  v 
coverage. 


1  The  National  Radio  Astronomy  Observatory  is  a  facility  of  the  National 
Science  Foundation  operated  under  cooperative  agreement  by  Associated 
Universities,  Inc. 

2  The  web  site  for  the  RRFID  is  located  at  http://rorf.usno.navy.mil/ 
rrfid.shtml  and  the  BVID  at  http://www.obs.u-bordeauxl.fr/BVID/ 


Observations  were  made  in  a  dual  frequency  bandwidth  synthesis 
mode  to  facilitate  delay  measurements  for  astrometry.  Observations 
in  this  mode  also  allow  imaging  at  both  frequency  bands.  Eight 
individual  frequency  bands  (IFs)  were  recorded  simultaneously, 
each  8  MHz  wide,  with  four  at  2.3  GHz  and  four  at  8.4  GHz  for  a 
total  of  32  MHz  in  each  frequency  band. 

The  data  were  correlated  with  the  VLBA  correlator  at  the  Array 
Operations  Center  in  Socorro,  New  Mexico.  The  correlated  data 
were  calibrated  and  corrected  for  time-  and  frequency-dependent 
phase  variations  using  NRAO’s  Astronomical  Imaging  Processing 
System  (AIPS;  Greisen  1998).  Initial  amplitude  calibration  for  each 
IF  was  accomplished  using  system  temperature  measurements  and 
gain  curves  generated  during  observations.  Fringe  fitting  was  done 
in  AIPS  using  a  solution  interval  equal  to  the  scan  duration  and  a 
point- source  model. 

Amplitude  calibration  was  improved  in  the  second  stage  by  using 
observations  of  sources  whose  core  flux  density  is  known  to  be 
>90  per  cent  of  the  total  flux  density.  For  a  precise  definition  of 
the  core,  refer  to  Section  3.1.  To  this  end,  a  single  amplitude  gain 
correction  factor  was  derived  for  each  antenna  based  on  fitting  a 
Gaussian  model  to  the  core  component.  Gain  correction  factors  were 
then  calculated  based  on  the  difference  between  the  observed  and 
the  model  visibilities.  Finally,  the  amplitude  gain  correction  factors 
were  applied  to  the  target  sources.  This  is  a  non-standard  procedure 
that  improves  the  overall  amplitude  calibration.  The  accuracy  of 
the  amplitude  calibration  determined  in  this  way  is  conservatively 
estimated  to  be  within  20  per  cent. 

Data  imaging  and  model  fitting  were  done  using  the  caltech  dif¬ 
ference  mapping  program  ‘difmap’  (Shepherd,  Pearson  &  Taylor 
1995),  after  inspecting  the  data  and  editing  out  obvious  bad  points. 
The  data  were  imaged  using  difmap  in  the  automatic  mode.  Gener¬ 
ally,  this  mode  fails  for  about  one-third  of  the  sources  which  have 
a  structure  that  is  complex  or  too  extended  for  the  automatic  script 
to  handle.  These  have  to  be  redone  by  hand  in  an  interactive  mode. 
For  the  southern  sources  in  this  sample,  almost  all  had  to  be  imaged 
by  hand  due  to  the  poor  u  —  v  coverage. 

difmap  combines  the  visibilities  in  all  four  2.3  GHz  receiver 
IFs  but  does  not  correct  for  spectral  index  effects.  It  was  assumed 
that  the  source  structure  variations  across  the  IFs  were  negligible 
(6  per  cent  variation  in  frequency  at  this  band).  Uniform  weight¬ 
ing  was  used  for  the  initial  phase  self-calibration  before  chang¬ 
ing  to  natural  weighting.  For  our  arrays,  uniform  weighting  gives 
more  weighting  to  the  longer  baselines,  whilst  natural  weighting 
gives  more  weighting  to  the  shorter  baselines.  Images  of  six  of  the 
31  sources  in  our  sample  are  shown  in  Fig.  2  with  the  remainder 
available  online  (see  supporting  information).  Images  of  additional 
epochs  can  be  found  in  the  form  of  contour  plots  at  the  RRFID  and 
BVID  websites. 

Generally,  circular  Gaussian  models  were  used  to  fit  the  u  —  v 
data  in  order  to  parametrize  the  source  morphology.  Like  imaging, 
model  fitting  is  an  iterative  process.  Elliptical  Gaussian  compo¬ 
nents  were  used  only  to  represent  the  core  component  or  a  very 
bright  jet  component  if  the  residuals  remaining  from  a  circular 
Gaussian  model  were  too  large  and  made  it  difficult  to  continue 
model  fitting  using  the  residual  map.  The  model  fits  generally  de¬ 
scribe  the  visibility  of  data  well,  but  these  models  may  not  be  unique 
because  of  incomplete  sampling  in  the  u  —  v  plane.  In  order  to  de¬ 
termine  the  suitability  of  a  source  as  a  calibrator,  we  then  determine 
the  amount  of  source  structure,  as  well  as  its  variation  with  time, 
using  several  different  methods  based  on  both  the  clean  compo¬ 
nents  from  the  difmap  imaging  and  on  the  parameters  of  the  fitted 
Gaussians. 
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Table  1.  Optical  and  radio  properties  of  the  sample.  The  positions  shown  above  are  the  most  recent  International  Celestial  Reference  Frame 
positions  (Fey,  Gordon  &  Jacobs  2009).  Score  is  the  flux  density  in  the  Gaussian  component  fitted  to  the  core  of  the  latest  epoch,  while  the  total 
flux  density  Stotai  is  the  sum  of  the  clean  components. 
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-03  50  04.6166450 

4.78 

4.73 

2008  Jan 

1908-201 

Quasar 

1.12* 

19  11  09.652  89198 

-20  06  55.108  9891 

2.10 

1.78 

2007  Dec 

1921-293 

Quasar 

0.35'' 

19  24  51.055  955  14 

-29  14  30.1210524 

6.18 

5.43 

2007  Dec 

1954-388 

Quasar 

0.63“ 

19  57  59.819  27470 

-38  45  06.355  7585 

2.75 

2.69 

2008  Jan 

1958-179 

Quasar 

0.65“ 

20  00  57.090444  85 

-17  48  57.6725440 

1.53 

1.56 

2007  Dec 

2052-474 

Quasar 

1.49“' 

20  56  16.359  818  74 

-47  14  47.627  6461 

1.37 

1.43 

2007  Dec 

2243-123 

Quasar 

0.63“ 

22  46  18.23197613 

-12  06  51.2774796 

2.57 

1.33 

2008  Jan 

2255-282 

Quasar 

0.93“ 

22  58  05.962  884  81 

-27  58  21.2567425 

0.95 

0.87 

2008  Jan 
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'  Jones  et  al.  (2009);  sdi  Serego -Alighieri  et  al.  (1994);  'Halpern,  Eracleous  &  Mattox  (2003);  “Browne,  Savage  &  Bolton  (1975). 


-90° 


Figure  1.  Sky  distribution  of  the  southern  sources  plotted  on  an  Aitoff 
equal-area  projection  of  the  celestial  sphere.  The  dotted  line  represents  the 
Galactic  plane,  while  the  dashed  line  is  the  ecliptic. 

3  ANALYSIS  AND  RESULTS 

Characterization  of  the  morphology  of  an  AGN  is  very  complex, 
and  no  single  metric  is  able  to  adequately  define  a  good  calibrator. 
Thus,  we  have  developed  a  number  of  complementary  approaches 


that  are  described  below.  A  weighted  combination  of  all  of  these 
metrics  is  eventually  used  to  classify  the  suitability  of  each  source 
to  be  a  high-resolution  phase  calibrator.  These  are  discussed  in 
the  following  sections. 

3.1  Core  flux  density 

A  good  calibrator  should  be  relatively  bright  at  the  frequency  of 
observation  to  be  easily  detectable.  It  should  also  be  stable  with 
minimal  flux  density  variation  over  time.  There  is  no  precise  def¬ 
inition  of  the  AGN  ‘core’  in  the  literature.  In  general,  the  bright, 
compact  flat  spectrum  feature  is  referred  to  as  the  core.  We  have 
confirmed  our  identification  of  the  cores  for  our  entire  sample  by 
establishing  that  these  structures  have  flat  or  inverted  spectra.  This 
was  done  by  using  the  8.4  GHz  data  that  are  observed  simultane¬ 
ously  with  the  2.3  GHz  observations  presented  here.  Here,  we  define 
a  core  flux  density,  Score,  for  each  source  and  for  each  epoch  as  the 
flux  density  of  the  Gaussian  component  fitted  to  the  core.  The  latest 
core  flux  density  as  defined  in  this  section  is  shown  for  each  source 
in  Table  1 .  From  the  core  flux  density,  we  computed  the  mean  core 
flux  density  ( S )  averaged  over  all  epochs  in  which  the  source  was 
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0727-115  at  2.309  GHz  2008  Jan  23 


I 

40  20  0  -20  -40 

Right  Ascension  (mas) 


1144-379  at  2.309  GHz  2007  Aug  01 


1908-201  at  2.292  GHz  2000  Oct  23 


1351-018  at  2.309  GHz  2006  Sep  13 


Figure  2.  Two  examples  each  of  class  ‘A’,  ‘B’  and  ‘C’  sources,  respectively.  Our  current  sample  contains  no  examples  of  classes  ‘D’  and  ‘E\ 
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Table  2.  Results  of  variability  analysis  in  the  flux  density,  core  fraction  and  radial  extent. 


Source  Epochs  Core  flux  density  Core  fraction  Weighted  radial  extent  Unweighted  radial  extent 


name 

S 

°core 

°core/  S 

C 

crc/C 

Rw 

°^W 

&rw/Rw 

^uw 

aR\jw 

crRim/R\jw 

0003-066 

24 

1.55 

0.31 

0.20 

0.64 

0.07 

0.11 

1.67 

0.40 

0.24 

6.17 

2.22 

0.36 

0104-408 

28 

1.26 

0.29 

0.23 

1.00 

0.04 

0.04 

0.93 

1.24 

1.33 

7.86 

17.13 

2.18 

0238-084 

16 

0.91 

0.20 

0.22 

0.53 

0.10 

0.19 

0.37 

0.10 

0.27 

17.34 

16.13 

0.93 

0336-019 

25 

1.92 

0.48 

0.25 

0.80 

0.08 

0.10 

1.00 

0.64 

0.64 

3.05 

2.32 

0.76 

0402-362 

17 

1.22 

0.11 

0.09 

1.00 

0.03 

0.03 

1.18 

0.26 

0.22 

8.11 

4.54 

0.56 

0454-234 

27 

1.73 

0.59 

0.34 

1.00 

0.02 

0.02 

0.46 

0.31 

0.67 

3.18 

1.78 

0.56 

0458-020 

26 

0.78 

0.31 

0.40 

0.68 

0.15 

0.22 

2.22 

1.02 

0.46 

16.30 

17.93 

1.10 

0727-115 

32 

2.37 

0.71 

0.30 

0.83 

0.05 

0.06 

1.15 

0.30 

0.26 

6.10 

1.77 

0.29 

0919-260 

17 

1.32 

0.83 

0.63 

0.73 

0.11 

0.15 

1.70 

0.34 

0.20 

6.46 

1.55 

0.24 

0920-397 

16 

1.06 

0.18 

0.17 

0.92 

0.11 

0.12 

1.81 

0.94 

0.52 

7.27 

2.98 

0.41 

1034-293 

27 

1.11 

0.30 

0.27 

0.88 

0.21 

0.24 

0.91 

0.68 

0.75 

5.83 

8.22 

1.41 

1124-186 

26 

0.90 

0.19 

0.21 

1.00 

0.02 

0.02 

0.25 

0.29 

1.18 

6.88 

16.50 

2.40 

1144-379 

23 

1.34 

0.39 

0.29 

0.83 

0.05 

0.06 

0.44 

1.79 

4.10 

9.61 

20.66 

2.15 

1145-071 

16 

0.69 

0.11 

0.16 

0.77 

0.10 

0.13 

1.33 

0.24 

0.18 

6.11 

2.26 

0.37 

1253-055 

3 

6.75 

0.27 

0.04 

0.65 

0.15 

0.23 

2.27 

0.25 

0.11 

12.5 

0.50 

0.04 

1255-316 

14 

1.20 

0.24 

0.20 

0.79 

0.11 

0.14 

2.68 

0.67 

0.25 

14.12 

11.86 

0.84 

1313-333 

17 

0.80 

0.28 

0.35 

0.71 

0.05 

0.07 

1.29 

0.40 

0.31 

7.82 

3.05 

0.39 

1334-127 

25 

2.12 

0.55 

0.26 

1.00 

0.04 

0.04 

0.84 

0.62 

0.74 

3.60 

1.87 

0.52 

1351-018 

13 

0.78 

0.25 

0.32 

1.00 

0.03 

0.03 

1.51 

2.93 

1.94 

2.55 

1.86 

0.73 

1424-418 

18 

1.69 

0.49 

0.29 

0.73 

0.08 

0.11 

3.42 

1.88 

0.55 

22.55 

11.07 

0.47 

1451-375 

14 

1.11 

0.30 

0.27 

0.88 

0.07 

0.08 

0.37 

3.71 

10.08 

9.31 

5.40 

0.58 

1514-241 

16 

1.76 

0.30 

0.17 

0.83 

0.05 

0.06 

3.39 

0.78 

0.23 

18.57 

3.90 

0.21 

1622-253 

24 

1.27 

0.38 

0.30 

0.83 

0.05 

0.06 

1.26 

0.72 

0.57 

6.35 

7.37 

1.16 

1741-038 

28 

3.18 

1.05 

0.33 

1.00 

0.03 

0.03 

0.33 

0.22 

0.66 

2.17 

1.24 

0.57 

1908-201 

23 

1.95 

0.45 

0.23 

0.77 

0.07 

0.09 

1.50 

0.45 

0.30 

5.60 

1.68 

0.30 

1921-293 

23 

7.42 

2.45 

0.33 

0.60 

0.09 

0.15 

2.30 

0.53 

0.23 

8.78 

8.17 

0.93 

1954-388 

21 

2.16 

0.54 

0.25 

1.00 

0.05 

0.05 

0.56 

0.52 

0.93 

3.52 

2.71 

0.77 

1958-179 

9 

1.16 

0.52 

0.45 

0.80 

0.04 

0.05 

0.17 

0.56 

3.29 

2.79 

2.29 

0.82 

2052-474 

10 

1.36 

0.34 

0.25 

0.81 

0.29 

0.36 

1.33 

2.09 

1.56 

10.45 

10.66 

1.02 

2243-123 

22 

1.50 

0.33 

0.22 

0.86 

0.06 

0.07 

2.00 

0.48 

0.24 

9.32 

1.77 

0.19 

2255-282 

20 

0.95 

1.11 

1.17 

0.77 

0.10 

0.13 

2.30 

2.05 

0.89 

6.58 

3.49 

0.53 

observed.  The  extent  to  which  the  core  flux  density  varies  over  time 
can  be  characterized  by  the  variability  index  (crcote/S)  of  the  core 
flux  density,  where  <7core  is  the  standard  deviation  of  the  core  flux 
density.  A  value  of  0.0  indicates  no  variation  over  time  (Table  2, 
Column  5). 

Fig.  3  shows  that  90  per  cent  of  the  sources  have  a  variability 
index  of  the  core  flux  density  below  0.4  and  23  per  cent  of  the 
sources  have  a  variability  index  below  0.2  (which  is  below  our 


Figure  3.  Distribution  of  the  core  flux  density  variability  index. 


estimated  calibration  uncertainties).  These  low  variability  indices 
indicate  that  most  of  the  core  fluxes  in  this  sample  are  stable. 

3.2  Core  fraction 

Following  Ojha  et  al.  (2004b),  we  define  a  core  fraction  as 

C  —  ‘Score  _cc/  ‘S'total?  (1) 

where  Scor&_cc  is  the  sum  of  the  flux  densities  of  clean  components 
within  one  synthesized  beam  of  the  brightest  pixel  and  Stotai  is  the 
sum  of  all  the  clean  component  flux  densities. 

It  provides  an  indication  of  how  point-like  a  source  is  and  also 
provides  a  way  to  track  source  structure  changes  from  epoch  to 
epoch.  The  average  core  fraction,  C,  was  computed  for  each  source 
over  all  epochs.  As  with  the  variability  index  of  the  core  flux  density, 
we  also  computed  the  variability  index  ( crc/C )  of  the  source  core 
fraction,  where  crc  is  the  standard  deviation  of  the  core  fraction. 
A  value  of  0.0  for  the  variability  index  of  the  source  core  fraction 
indicates  no  variation  over  time.  The  results  are  shown  in  Table  2, 
Column  8. 

The  distribution  of  the  variability  index  of  the  core  fraction  is 
shown  in  Fig.  4,  while  the  distribution  of  the  core  fraction  is  shown 
in  Fig.  5.  The  mean  core  fraction  for  all  the  sources  is  83  per  cent 
with  a  standard  deviation  of  12  per  cent.  All  sources  have  a  mean 
core  fraction  with  the  variability  index  below  0.4.  27  sources  have  a 
variability  index  between  0  and  0.2  and  the  remaining  four  between 
0.21  and  0.4.  In  general,  the  southern  sample  sources  are  very 
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Figure  4.  Distribution  of  the  variability  index  of  the  core  fraction.  Note  that 
all  sources  have  a  variablity  index  <0.4. 
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Figure  5.  Distribution  of  the  core  fraction.  This  measure  was  used  instead 
of  the  variability  index  as  it  shows  slightly  higher  variation  than  the  core 
fraction  variability  index. 

compact  with  little  variation.  For  this  reason,  we  used  the  actual 
core  fraction  and  not  the  variability  index  to  classify  the  sources. 


3.3  Flux-weighted  radial  extent 


The  previous  two  metrics  are  primarily  measures  of  core¬ 
dominance.  This  and  the  following  metric  provide  complementary 
information  by  quantifying  the  radial  extent  of  the  sources  in  the 
Southern  sample. 

The  flux-density-weighted  radial  extent  (Ojha  et  al.  2004b)  is 
defined  as 


Rw  — 


Ei  sin 

E;  si  ’ 


(2) 


where  Rw  is  in  units  of  milliarc seconds  and  rt  is  the  radius  at  which 
the  zth  clean  component  has  the  flux  density  St. 

The  mean  (Rw)  and  standard  deviation  were  also  calculated  as 
was  the  variability  index  (crRw/Rw),  where  <jRw  is  the  standard 
deviation  of  the  weighted  radial  extent,  see  Table  2,  Column  1 1 . 

The  variability  index  of  the  weighted  radial  extent  is  widely 
distributed  as  shown  in  Fig.  6,  with  10  per  cent  of  the  sources 
having  an  index  between  0  and  0.2,  showing  minimal  variation  in 
the  weighted  radial  extent.  32  per  cent  of  the  sources  have  an  index 
between  0.21  and  0.40,  while  13  per  cent  have  an  index  between 


Figure  6.  Distribution  of  the  variability  index  of  the  weighted  radial  extent. 
Seven  sources  with  the  variability  index  greater  than  1  are  not  shown  in  the 
histogram. 

0.41  and  0.60.  16  per  cent  of  the  sources  have  an  index  between 
0.61  and  0.80,  while  the  rest  (29  per  cent)  have  an  index  greater 
than  0.81.  This  shows  a  high  degree  of  variability  in  the  weighted 
radial  extent. 


3.4  Unweighted  radial  extent 

This  is  the  radial  extent  within  which  95  per  cent  of  the  flux  density 
of  the  source  is  contained.  This  provides  a  measure  of  how  ex¬ 
tended  a  source  is.  The  mean  (/?uw)  and  standard  deviations  were 
calculated  as  was  the  variability  index  (crRim/ R\jw),  where  crRim  is 
the  standard  deviation  of  the  unweighted  radial  extent  and  /?uw  is 
in  units  of  milliarc  seconds  (Table  2,  Column  14). 

Like  the  variability  index  of  the  weighted  radial  extent  (see 
above),  the  variability  index  of  the  unweighted  radial  extent  is  also 
widely  distributed.  Fig.  7  shows  that  6  per  cent  of  the  sources  have 
an  index  below  0.2,  23  per  cent  have  an  index  between  0.21  and 
0.40,  while  29  per  cent  have  an  index  between  0.41  and  0.60.  About 
10  per  cent  have  an  index  between  0.61  and  0.80,  while  35  per  cent 
have  an  index  greater  than  0.80,  showing  very  high  variability  in  the 
unweighted  radial  extent.  As  discussed  below,  this  high  variability 
(in  both  the  weighted  and  unweighted  radial  extents)  is  largely  a 


Variability  Index 


Figure  7.  Distribution  of  the  variability  index  of  the  unweighted  radial 
extent.  Six  sources  with  the  variability  index  greater  than  1  are  not  shown 
in  the  histogram. 
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result  of  low  brightness  ‘jet’  features  that  are  detected  in  some 
epochs  but  not  others,  usually  due  to  differing  sensitivities  of  the 
observing  array  rather  than  a  real  variation  in  the  source. 

4  SOURCE  CLASSIFICATION 

To  capture  the  complexity  of  determining  whether  a  radio  source 
is  a  suitable  phase  calibrator  for  a  radio  interferometer,  we  have 
developed  the  set  of  complementary  metrics  described  above.  Here, 
we  discuss  how  we  combine  the  information  these  metrics  provide 
to  arrive  at  a  classification  scheme  for  potential  calibrators.  We  ex¬ 
plain  our  reasons  for  weighting  the  different  metrics  as  we  have 
done  and  present  our  recommendations  based  on  this  weighting 
scheme.  Since  such  a  scheme  is  necessarily  somewhat  subjective,  it 
should  be  considered  as  reasonable  and  useful  rather  than  a  defini¬ 
tive  classification. 

The  four  metrics  we  used  to  classify  the  calibrator  source  quality 
were  core  fraction,  core  flux  density  variability  index,  weighted  and 
unweighted  radial  extents.  Each  of  these  metrics  is  first  assigned 
a  score,  with  the  scores  as  shown  in  Tables  3  and  4.  The  scores 
for  the  individual  metrics  were  chosen  so  as  to  sum  to  100  for  a 
perfect  calibrator  source.  We  also  chose  to  score  the  core  fraction 
and  core  flux  density  variability  so  as  to  give  them  higher  weight, 
each  having  a  maximum  score  of  35,  while  the  two  radial  extents 
have  the  maximum  score  of  15.  We  use  a  lower  maximum  score 
for  the  radial  extents  as  they  are  highly  sensitive  to  small  epoch- 
to-epoch  variations  in  the  signal-to-noise  ratio  of  the  images.  Faint 


Table  3.  Score  distributions  of  the  variability  indices  of 
the  source  core  fraction  and  core  flux  density. 


Core  flux  density  Score 
variability  index 
(per  cent) 

Core  fraction 
variability  index 

Score 

>1 

0.00 

0.91-1.00 

3.5 

91-100 

35.0 

0.81-0.90 

7.0 

81-90 

31.5 

0.71-0.80 

10.5 

71-80 

28.0 

0.61-0.70 

14.0 

61-70 

24.5 

0.51-0.60 

17.5 

51-60 

21.0 

0.41-0.50 

21.0 

41-50 

17.5 

0.31-0.40 

24.5 

31-40 

14.0 

0.21-0.30 

28.0 

21-30 

10.5 

0.11-0.20 

31.5 

11-20 

7.0 

0.01-0.10 

35 

1-10 

3.5 

Table  4.  Distribution  of  scores  for  the 
variability  index  of  the  weighted  and 
unweighted  radial  extents. 


Index  Score 


>1  0.00 

0.91-1.00  1.5 

0.81-0.90  3.0 

0.71-0.80  4.5 

0.61-0.70  6.0 

0.51-0.60  7.5 

0.41-0.50  9.0 

0.31-0.40  10.5 

0.21-0.30  12.0 

0.11-0.20  13.5 

0.01-0.10  15 


extended  features  may  be  detected  only  at  some  epochs  depend¬ 
ing  on  the  signal-to-noise  ratio,  resulting  in  a  large  variation  in  the 
radial-extent  measures.  Such  low  surface  brightness  features  gen¬ 
erally  have  negligible  impact  on  the  usefulness  of  the  source  as  a 
calibrator.  In  any  case,  for  a  lower  resolution  array  like  MeerKAT 
these  extended  structures  will  be  embedded  in  the  main  central 
component. 

Based  on  the  overall  score  summed  up  from  this  weighting 
scheme  each  source  falls  into  one  of  the  five  classes  ‘A’  through  ‘E’ . 
A  source  falls  into  class  A’  if  its  overall  score  is  between  80  and 
100  per  cent,  class  ‘B’  if  the  score  is  between  60  and  80  per  cent, 
class  ‘C’  if  it  scores  between  40  and  60  per  cent,  class  ‘D’  for  any 
score  between  20  and  40  per  cent  and  class  ‘E’  for  sources  scoring 
between  0  and  20  per  cent.  We  propose  the  following  classifications. 

A  -  excellent  calibrator  (score  between  80  and  100). 

B  -  very  good  calibrator  (score  between  60  and  80). 

C  -  good  calibrator  (score  between  40  and  60). 

D  -  use  with  caution  (score  between  20  and  40). 

E  -  unsuitable  as  a  calibrator  (score  below  20). 

In  our  sample  (Table  5),  we  have  nine  class  A’  sources,  19  class 
‘B’  sources  and  three  class  ‘C’  sources.  There  are  no  sources  in 
classes  ‘D’  and  ‘E’.  Thus,  all  the  sources  in  our  sample  would  be 
‘good’  calibrators  and  all  but  three  sources  are  likely  to  be  ‘very 
good’  calibrators. 


Table  5.  Classification  of  the  sources.  Wi:  core  flux  density;  W2: 
core  fraction;  W3:  unweighted  radial  extent;  W4:  weighted  radial 
extent. 


Source 

name 

Wi 

Score 

W2  w3 

W4 

Total  score 
(per  cent) 

Class 

0003-066 

31.5 

31.5 

10.5 

12.0 

85.5 

A 

0104-408 

28.0 

35.0 

0.0 

0.0 

63.0 

B 

0238-084 

28.0 

31.5 

1.5 

12.0 

73.0 

B 

0336-019 

28.0 

35.0 

4.5 

6.0 

73.5 

B 

0402-362 

35.0 

35.0 

7.5 

12.0 

89.5 

A 

0454-234 

24.5 

35.0 

7.5 

6.0 

73.0 

B 

0458-020 

24.5 

31.5 

0.0 

9.0 

65.0 

B 

0727-115 

28.0 

35.0 

12.0 

12.0 

87.0 

A 

0919-260 

14.0 

31.5 

12.0 

13.5 

71.0 

B 

0920-397 

31.5 

31.5 

9.0 

7.5 

79.5 

B 

1034-293 

28.0 

28.0 

0.0 

4.5 

60.5 

B 

1124-186 

28.0 

35.0 

0.0 

0.0 

63.0 

B 

1144-379 

28.0 

35.0 

0.0 

0.0 

63.0 

B 

1145-071 

31.5 

31.5 

10.5 

13.5 

87.0 

A 

1253-055 

35.0 

28.0 

15.0 

13.5 

91.5 

A 

1255-316 

31.5 

31.5 

3.0 

12.0 

78.0 

B 

1313-333 

24.5 

35.0 

10.5 

10.5 

80.5 

A 

1334-127 

28.0 

35.0 

7.5 

4.5 

75.0 

B 

1351-018 

24.5 

35.0 

4.5 

0.0 

64.0 

B 

1424-418 

28.0 

31.5 

9.0 

7.5 

76.0 

B 

1451-375 

28.0 

35.0 

7.5 

0.0 

70.5 

B 

1514-241 

31.5 

35.0 

12.0 

12.0 

90.5 

A 

1622-253 

28.0 

35.0 

0.0 

7.5 

70.5 

B 

1741-038 

24.5 

35.0 

7.5 

6.0 

73.0 

B 

1908-201 

28.0 

35.0 

12.0 

12.0 

87.0 

A 

1921-293 

24.5 

31.5 

1.5 

12.0 

69.5 

B 

1954-388 

28.0 

35.0 

4.5 

1.5 

69.0 

B 

1958-179 

21.0 

35.0 

3.0 

0.0 

59.0 

C 

2052-474 

28.0 

24.5 

0.0 

0.0 

52.5 

C 

2243-123 

28.0 

35.0 

13.5 

12.0 

88.5 

A 

2255-282 

00.0 

31.5 

7.5 

3.0 

42.0 

C 
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5  CONCLUSIONS 

We  model  fitted  up  to  32  epochs  of  observations  (average  of  20 
epochs)  for  each  of  the  31  sources  in  our  sample  which  was  se¬ 
lected  from  the  RRFID  kinematic  survey  (Piner  et  al.  2007)  and  de¬ 
termined  their  suitability  as  phase  calibrators.  While  the  kinematic 
survey  looks  more  at  proper  motion  in  the  sources  at  8.4  GHz,  this 
paper  concentrates  on  the  morphological  properties  of  the  sources  at 
2.3  GHz,  a  frequency  more  relevant  to  emerging  Southern  hemi¬ 
sphere  arrays  like  MeerKAT  and  ASKAP. 

We  have  developed  a  method  to  classify  radio  sources  according 
to  their  suitability  as  phase  calibrators  for  radio  interferometers.  We 
first  characterize  a  source  by  calculating  several  metrics  which  give 
measures  of  the  degree  to  which  the  core  dominates  the  source, 
and  the  degree  of  variability,  both  in  flux  density  and  the  degree  to 
which  the  source  is  extended.  These  metrics  are  then  combined  to 
give  the  source  a  total  score,  which  is  used  to  assign  the  source 
to  one  of  the  five  classes  of  suitability  as  a  calibrator.  All  31  sources 
in  our  sample  were  classified  as  ‘good’  calibrators  with  28  classified 
as  ‘very  good’  or  better. 
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SUPPORTING  INFORMATION 

Additional  Supporting  Information  may  be  found  in  the  online  ver¬ 
sion  of  this  article: 

Figure  2.  Images  of  the  3 1  sources  in  our  sample. 

Please  note:  Wiley-Blackwell  are  not  responsible  for  the  content  or 
functionality  of  any  supporting  materials  supplied  by  the  authors. 
Any  queries  (other  than  missing  material)  should  be  directed  to  the 
corresponding  author  for  the  article. 
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